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Abstract 



The dynamical symmetry breaking in a two-field model is studied by nu- 
merically solving the coupled effective field equations. These are dissipative 
equations of motion that can exhibit strong chaotic dynamics. By choos- 
ing very general model parameters leading to symmetry breaking along one 
of the field directions, the symmetry broken vacua make the role of transi- 
tory strange attractors and the field trajectories in phase space are strongly 
chaotic. Chaos is quantified by means of the determination of the fractal di- 
mension, which gives an invariant measure for chaotic behavior. Discussions 
concerning chaos and dissipation in the model and possible applications to 
related problems are given. 
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I. INTRODUCTION 



The study and understanding of the dynamics of fields are a timely subject and of 
broad interest, with applications in diverse areas like in particle physics, cosmology and in 
condensed matter (for a recent review, see and references therein). Additional interest on 
the subject comes from the fact that many of the theoretical ideas and models can be tested 
in ongoing experiments, as those been performed in condensed matter systems, and in the 
future ones, in the entering of operation of the RHIC and LHC heavy ion colliders, which will 
be able to probe possible new phenomena at the QCD scale and space based experiments, 
which will be putting on test different cosmological models. It is then becoming urgent the 
detailed investigation of the underline field dynamics that may be common to all these very 
different areas of physical research. 

In this paper we are mainly concerned with the connection between the development of 
strong nonlinearities in the time evolving system of equations of motion of a given field theory 
model and the possible chaotic behavior associated to them. Lets recall that in symmetry 
breaking phase transitions we are usually interested in the study of the evolution of a given 
order parameter, for example the magnetization in spin systems in statistical physics or a 
vacuum expectation value of some scalar field (the Higgs) in particle physics, which gives 
a measure of the degree of organization of the system at the macroscopic level. However, 
at the microscopic level disorder is related to the nonlinearities and fluctuations responsible 
to chaotic behavior in the system and these chaotic motion phenomena can reflect in a 
nontrivial time dependence of the macroscopic quantities and, therefore, influencing all the 
dynamics of the system. This is clear once several properties of the system at longer times 
are closely related to the microscopic physics, like relaxation to equilibrium, phase ordering, 
thermalization and so on. Thus we expect that chaos will be not only an important ingredient 
in determining the final states of a given system but also in how it gets there. 

Previous studies on chaos in field theory have mostly emphasized chaotic behavior in 
gauge theory models (for a review and additional references and applications, see 0). In 
special, in homogeneous Yang-Mills-Higgs models we can reduce the system of classical 
equations of motion to ones analogous to those of nonhnear coupled oscillators, which is 
well known to exhibit chaotic motion. 

While in all the previous works on chaotic dynamics of fields dealt with (conservative) 
Hamiltonian systems, here we will be mainly concerned with the effective field evolution 
equations, which are known to be intrinsically dissipative [0-^] and, therefore, the dynamical 
system we will be studying is non-Hamiltonian. 

Thus, we will be concerned with the influence of field dissipation, due to field decaying 
modes, on the degree of chaoticity of the field dynamics. Typically we expect that dissipation 
damps the fluctuations on the system and consequently tends to suppress possible chaotic 
motions and makes field trajectories in phase space to tend faster to the system asymptotic 
states. On the other hand there are well known examples of dynamical systems, as the 
Lorenz system [|T^], which are dissipative ones and at the same time they display a very rich 
evolution on phase space in which, under appropriate system parameters, field trajectories 
may be lead towards strange attractors. The verification of the same properties in a model 
motivated by particle physics would be a novel result with possible consequences to, for 
example, particle physics phenomenology and cosmology. 
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We study chaos in our dynamical system of equations by means of the measure of the 
fractal dimension (or dimension information) [for a review and definitions, see e.g., [[10|] ], 
which gives a topological measure of chaos for different space-time settings and it is a 
quantity invariant under coordinate transformations, providing then an unambiguous signal 
for chaos |Tl|,|T2|. The method we apply in this work for quantifying chaos will then be 
particularly useful in our planned future applications of our model and it extensions to a 
cosmological context, in which case other methods may be ambiguous, like, for example, the 
determination of Lyapunov exponents, which does not give a coordinate invariant measure 
for chaos, as discussed in |TT],0. Also, other methods for studying chaotic systems, like for 



example by Poincare sections, are not suitable in the case we are interested here, in which 
chaos is a transitory phenomenon as we will see. 

This work is organized as follows. In Sec. II we introduce the model and discuss its 
general properties at the classical level. In Sec. Ill we obtain the one-loop effective equations 
of motion for the fields and we determine the general form of the nonlocal (non-Markovian) 
dissipative kernels. In Sec. IV we then discuss the vahdity of the one-loop approximation and 
the Markovian approximation for the dissipative kernels appearing in the one-loop effective 
equations. We couple our system of fields to a set of N other fields making up the bath 
(or environment), in which the system evolves, and the large N behavior of the various 
important quantities is determined. In Sec. V we then present the Markovian form for the 
effective equations of motion and our main numerical results, where we determine the fractal 
dimension. In Sec. VI our concluding remarks are given. 



II. THE TWO INTERACTING SCALAR FIELDS MODEL 

The model we will study consists of two scalar fields in interaction with Lagrangian 
density given by 

All coupling constants are positive and > 0, but we choose < 0, such that we allow 
for spontaneous symmetry breaking in the \l/-field direction. Additionally, note that from 
the above Lagrangian, that for values of $ larger than a where = \m'^\/g'^, there 
is no symmetry breaking in the \l/-field direction. Thus, for example, if we have an initial 
state prepared at $ > ^cr, the \1/ field will move towards zero and remain around that 
state till eventually $ crosses below the critical value inducing a (dynamical) symmetry 
breaking in the \l/-field direction, after which the fields evolve to their vacuum values at 
(^)^ = ±V6\m^\/^ and = 0. 

The classical equations of motion for field configurations ($) = 0c and (\I') = ipc can be 
readily be obtained from Eq. (^.11): 

□0, + mj0, + ^(Pl + ^?Vc^e = , (2.2) 



2 4! 
2 4! 
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(2.1) 
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+ mj^, + ^iPl + = . (2.3) 



For homogeneous fields <^c and ipc (in wliicli case Eqs. ( p.2| ) and ( |2.3| ) are equivalent to 
the equations of motion of two particles with quartic potentials and quadratic interaction 
between them) the above equations are well known to lead to chaotic trajectories in phase 
space. The chaotic behavior of very similar classical equations have been studied recently 



in Ref. [ITl] for a model with Z2 x Z2 symmetry and the corresponding dynamical system 



shown to be chaotic for symmetry breaking in one of the field directions. 

III. THE ONE-LOOP EQUATIONS OF MOTION 

In the following we study how quantum effects, which will be responsible for the appear- 
ance of dissipative dynamics in the system evolution, will influence this dynamical symmetry 
breaking process and we identify and quantify possible chaotic behavior in the system evo- 
lution, function of the fields dissipations. 

Roughly speaking, chaos means extreme sensitivity to small changes in the initial con- 
ditions. Due to nonlinearity, fluctuations in the initial conditions of chaotic systems evolve 
such that they can completely alter the asymptotic outcome of the unperturbed trajectories 
in phase space. We here quantify the chaoticity of the system by measuring the fractal 
dimension. The fractal dimension is associated with the possible different exit modes under 
small changes of initial conditions and it gives a measure of the degree of chaos of a dynami- 



cal system |T0|. The exit modes we refer to above are one of the symmetry breaking minima 
in the ^'-field direction, which are attractors of field trajectories in phase space. The method 
we employ to determine the fractal dimension is the box-counting method, whose definition 
and the specific numerical implementation we use here have been described in details in Ref. 

Quantum corrections are taken into account in the effective equations of motions (EOMs) 
by use of the tadpole method |TB[. Let 0c and ifjc be the expectation values for $ and 



respectively. Splitting the fields in ( p.l|) in the expectation values and fluctuations, 

$-^0c + and *^Vc + ^, (3.1) 

where ($) = 0c and (\E') = ipci the EOMs for 0c and ipc are obtained by imposing that 
(0) = and {ip) = 0, which lead to the condition that the sum of all tadpole terms for 
each field vanishes. Restricting our analysis of the EOMs to homogeneous fields (0c = 0c(^), 
ipc = '0c(t)), thus, at the one-loop order we can write the following EOMs for 0c and tpc, 
respectively, as[| (where overdots mean time derivatives) 

0. + m50c+^0^+^7Vc^c' + ^0c(0^) 

+ /0c(^2) + 2/V'c(0^) = O (3.2) 



^Note that the mixed field averages in ( |3.2| ) and (3.3) can only be treated within perturbation 
theory. 
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and 



,2„;, /i2\ I o„2. 



(3.3) 



where (0^) and (■?/'^) are given in terms of the coincidence hmit of the (causal) two-point 
Green's functions G^^{x,x') and G^'^{x,x'), which are obtained from the (1, l)-component 
of the real time matrix of full propagators which satisfy the appropriate Schwinger- Dyson 
equations (see, e.g., Refs. [H and for further details): 



A. 



G^{x,x') 



+ / d'^zTj^{x, z)G(j,{z, x') = i6{x,x') 



(3.4) 



and 



GJx,x') 



+ J d'^z'E^{x, z)G^{z, x) =i6{x,x') , (3.5) 

where S(^(x,x') and S the self-energies for $ e \E', respectively. The momentum- 

space Fourier transform of G{x,x') (for both fields) can be expressed in the form 



G{x, x') 



, , ^ .q.(x-x') / G++(q, t - t') G+-(q, t - t') 
' {2Tif \G-^{c^,t-t') G-{c^,t-t') 



(3.6) 



where 



G'++(q, t-t') = G>(q, t - t')9{t - t') + G<(q, t - t')e{t' - t) 
G— (q, t-t') = G> (q, t - t')e{t' -t) + G< (q, t - t')9{t - t') 
G+-(q,t-t') = G<(q,t-t') 
G-+(q,t-t') = G>(q,t-t') , 



(3.7) 



and the fully dressed two-point functions, at a finite temperature T = are given byQ 
G>(q,t - t') = ^ {[1 + n(cu - ^r)]e-'("-'^)(*-*') + n(cu + tT)e'^'^+'^^^'-''^} e{t - t') 



+ 



^[[l + n{uj + ir)]e-^('"+^r)(*-*') + n{uj - iT)e'^^~'^^^'-''^] e{t' - t) 



G<{q,t-t')=G>{ci,t'~t) 



(3.8) 



where n{u) is the Bose distribution, u = u^q) is the particle's dispersion relation and T is 
the decay width, defined as usual in terms of the field self-energy by 



^R.O.R. thanks I. Lawrie for pointing him the correct form for these expressions. 
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ImE(q, u) 

2ij 



(3.9) 



Assuming the couplings g,X(f„X^ <^ 1, so that perturbation theory can be consistently 
formulated and subleading terms can be neglected, then by perturbatively expanding the 
field averages in Eqs. (|3.2| ) and (|3.3|) , we can write for the EOMs 



+2g^(t)c{t) 



dt' 



dt' 



(27r)3 



(2vr) 



Im 
Im 



Gl+(q,t-0 



+8^7Vc(t) / dt'<P,{t')^lj,{t') 



d^q 
(2^ 



Im 



and 



dt' 



dt' 



d^q 
d^q 



Im 
Im 



G^+{q,t-t') 
G^+{q,t-t') 



+8g%{t) / dt'Mt')Mt') 



d^q 

(2^ 



Im 



(27r)3 

G;+{q,t-t')G;+{q,t-t' 



where 



A. 



^2 , A^,^,2, 



(3.10) 



(3.11) 



(3.12) 



with (. . .)o meaning fields independent averages. The nonlocal terms of the type appearing 
in the above equations have been shown in Refs. to lead to dissipative dynamics in the 
EOMs. This can be made more explicit by an appropriate integration by parts in the time 
integrals in Eqs. ( p.lO|) and (|3.11| ) (see Ref. |§|) to obtain the result 



A 



+Ut) / dt'Mt')(P,it')Fiit, t') + 0,(t) / dt'ij,{t')i:,{t')F^{t, t') 

+z^,(t) /* dt' k(t')4(t') + Mt')Mt')] F^it, t') = 0, (3.13) 



and 
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+<Pc{t) / dt' Mt')Mt') + Mt')Mt') F^{t,t') = 

where the dissipative kernels Fi, F2, F3 and F4 are given by 

d^q 



(3.14) 



Fi{t,t') 
F2M 
F,M 
Fi{t,t') 



dt' lS^\ A?, Tm 



3 1 V 



(2vr)= 



G++(q,t-t')J +Vlm[G++(q,t-t': 
G++(q, t - t')] ' + AgHm [G++(q, t - t' 
G++(q, t - t')] ' + Av,Im [G++(q, t - t' 
G++(q,t-tOG^+(q,t-t'^ 



(3.15) 



and A^, A^ and g'^ (the one- loop effective coupling constants) in Eqs. (|3.13| ) and (|3.14|) are 
given by 



A. = A,+ /rft'/^l^Im 



A^ — A^ + / dt' 



(27r)3 [ 2 



G++(q,t-t')l +2(7^Im[G++(q,t-t' 



G,r(q,t-t')l +2(7^Im[G++(q,t-t' 



t'=t 



f = g' + g'ldt' ^ 



A^Im 



(27r)3 

+ 8/lm [G++(q, t - t')G;+(q, t - t')] } 



G:+(q,t-t') 



+ A^Im 



Gl+(q,t-t') 



(3.16) 



IV. COUPLING TO BATH DEGREES OF FREEDOM AND AN APPROXIMATE 

FORM FOR THE EOMS 



The two coupled nonlocal EOMs, Eqs. ( p.l3|) and (|3.14|) (or equivalently, Eqs. (p.lO|) 
and ( p.ll|) ) are too complicated to directly numerically work with them. The main difficulty 
in handling these equations comes from the dissipative like terms in Eqs. (|3.13| ) and ( p. 141) , 
which have the non-Markovian kernels shown in Eq. ( p.l5| ). If we suppose that there is a 
Markovian limit for those kernels, then by using Eqs. (|3.7|) and (|3.8| ) in Eq. (|3.15| ), we find 
that at T = the kernels diverge logarithmic (a result also found in Ref. @] for the single 
scalar field case). However, in the high temperature limit (F/a; <C 1 and F/T -C 1) it has 
been argued in Ref. that a Markovian approximation exists and a finite result for the 
dissipation coefficients can be found. Such an approximation for the kernels given in Eq. 
( p.l5| ), we can also find here for our specific problem, in which case we may then, in principle, 
write for the dissipation terms in Eqs. ( p.l3| ) and ( |3.14| ) the approximate expressions 
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dt' 



dt' 



Mt) / dt'Mt')Mt')Fi{t,t') 
Mt) f dt'Mt')Mt')F2{t,t') 

J —oo 

0e(t) r dt'Mt')Mt')F3{t,t') 

J —oo 

Mt) f dt'Mt')Mt')F3{t,t') 

J —oo 

Mt')Mt') + Mt')Mt')] F4{t,t') 



(i)c{t)ipcit)(pcit)r]3 , 

Mt) [<l>c{t)Mt) + Mt)Mt) 



V4 
V4 



(4.1) 



where rii,'r]2,ri3 and 774 in the above expressions are given by (using Eqs. ( p.7| ) and ( pl8 



m 



Eqs. ( p.l3| ) and ( p.l4| ), and in the high temperature approximation T^/T,T^/T -C 1, with 
T^/u^, T^/uj^, < 1 and 7^ m^) 



112 ~ — P 



d^q n^{l + n^) 



PI 



(2vr) 



d^q n^{l + n^) 
(2vr)3 ulT^ 



d^p n^{l + n^) ^ ^X^^ f d^p n^{l + n^) 



r/4 ~ 8g^f3 J 



(27r)3 
d^q 



1 



(2vr) 



(27r)3 (c.2 - u;iy 



/pa p2 \ 

K(l + n^)T^ + n^(l + n^)T^] + O (7^4, (7^4 . (4.2) 



Similar time non-locahties as the ones appearing in Eqs. (|3.1CI|) and (|3.11|) have also 
been dealt with in Refs. by using an adiabatic (or sudden) approximation for the fields. 
The dissipation coefficients obtained by that approximation are just the same as the ones 
given in ( |4.2| ). As shown in Ref. the adiabatic approximation for the nonlocal kernels is 
a consistent approximation in the case the fields are in an overdamped regime. The strong 
field dissipation responsible for the overdamped regime can be attained by coupling both $ 
and \E' fields to a large set of other fields making the bath in which $ and \l/ evolve, which 
then enlarges the number of field decay channels available for both fields. This idea was 
used in Ref. for the construction of an alternative infiationary model. 



A. The Validity of the One-Loop and Markovian Approximations 

The consideration of the coupling of $ and \E' to a large set of bath field degrees of 
freedom is particularly relevant to the validity of the Markovian approximation taken for 
the dissipative kernels and in order to make a clear assessment of the regime of validity 
not only of this approximation but also for the validity of the one-loop approximation we 
used to derive the EOMs. The study of the validity of these approximations are clearly 
important to our study which is the study of chaotic behavior in the dynamics of our system 
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of equations, and by making clear that chaos in our dynamical system is not just an artifact 
of the approximations taken. In order to address these important issues, let us consider the 
coupling of both $ and \1/ to a set of (scalar) fields x making the bath, which are coupled 
in the following way: 



N N 



flY.^\l + flY.'^'xl ■ (4.3) 



i=l i=l 



We next study how our main quantities scale with the various coupling constants and 
in the large N limit]^. This analysis will allow us to make a qualitative and clear assessment 
of our main approximations. 

Let us first consider that the various coupling constants we have scale with in the 
following way: 



~ ~ / ~ — = , (4.4) 



and 

with a <^ 1. The choice for the coupling constants taken above is consistent in the sense that 
higher order quantum corrections to quantities like the field masses and coupling constants 
will not blow up in the large N limit. For instance, as — oo the effective masses m| and 
scale as (in the high temperature limit) 

= mj + 0{\^T') + 0{g'T') + 0{N flT^) ^ + 0{a'T') , 
= mj + 0{\^T') + 0{g'T') + 0{N f^T^) ^ mj + 0{a'T^) , (4.6) 

and the effective coupling constants go like 

„2 



A,^ = A^ + 0{Xl) + 0{g') + 0{Nfl) ^X^ + O 



A^ = A^ + 0(AJ) + 0{g') + 0{Nfl) ^ + O 



a 
N 



f=g' + 0{X^g') + 0{X^g') + 0{g') + 0(Ar/|/|) ^g^ + 0[j^j . (4.7) 

We next determine the scaling of the fields decay widths V^, and for those of the 
bath, r^-, which determines the time-scale for collisions for the system in interaction with 
the (thermal) bath. We begin by noticing that the imaginary part of the field self-energies 



^R.O.R. deeply thanks S. Jeon for correspondence regarding the large N contributions to dissipa- 
tion and the consistency of the adiabatic approximation in the large N limit. 
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contributing to the decay widths are typically dominated by momenta |q| ~ T [ll5|,|6|. We 
can then find for the decay widths the relations: 



a2 



= 0{\lT) + 0{g'T) + 0{NflT) ^ O 



N 



- Oif^T) + 0{flT) ^ O ( ) . (4.8) 



Using these in the expressions for the dissipative kernels, Eq. ( p.l5| ), we can then show that 
the dissipation coefficients 771, ... ,774, appearing in Eq. ([4.1D, have magnitude given by 



r/3 




Therefore, the above estimates show that, for <^ 1 and in the large N limit, the dis- 
sipation kernel F4 and its related dissipation coefficient in the Markov approximation are 
negligible compared with the others one. The behavior of the remaining coefficients repro- 
duce the typical behavior found previously for the dissipation coefficients in Refs. and . 
Additionally, we can also find that higher loop contributions to dissipation (for instance the 
ones coming from two-loop diagrams and higher) are all at most of order (^(A^r,^), O(A^r^), 
0{g'^T^ + g^T^) or 0{Nf^T^), which are all subleading in the perturbative regime. All these 
simple qualitative estimative allow us then to assess the validity of both the Markov and 
one- loop approximations used in this study. 



V. THE DYNAMICAL SYSTEM AND CHAOTIC BEHAVIOR 



From the results obtained in the last Section, we can then write Eqs. ( |3.13|) and (|3.14|) 



in the local form shown below and more suitable for the numerical analysis as: 

k + ml(pc + ^(Pl + f(i)c^l + Vi(p% + V3(pc^cipc = (5.1) 



and 
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where ~ A,^, ~ A^ and ~ g'^ denote the renormahzed couphngs ^ Viy V2 and 773 
denote the dissipation coefficients and we have neglect the term involving 774 due to the 
estimates shown in ( ^.91 ) . Although these coefficients can be explicitly evaluated, however, 
we here refrain ourselves from an explicit evaluation of these terms, which, as shown in the 
previous Section, depend on the various parameters of the model, temperature and on the 
number of other field degrees making the environment that $ and may be coupled to. 
In fact, the magnitude of the dissipation terms may be controlled by these additional field 
couplings, as shown in the simple estimates obtained in the previous Section. Therefore, for 
the sake of simplicity we just take rji, rj2 and 773 as additional free constant parameters. 
Next, let us define the following constants 



■2 



ml 
6|m^| 

-9' 



a^A^ 

in term of which we can define the dimensionless variables: 



A. = , (5.3) 



— a2 



X — \l Aw,^— <4ic 1 



w = —^= 'j/'c , (5.4) 



and by also rescaling time and dissipation coefficients as t' = y 6|m^| t, rji = rj^ a^y A^/ {'^)v, 
V2 = Vy y^tt'/ {^)v and ?73 = rjxy a^JX^/ {^)vi respectively, we can then write Eqs. (|0|) and 



^2^ in terms of the following dimensionless ffist-order differential system of equations: 



X = z 
y = w 

3 

w = -- — - Gxy-r]yyw- r]xyxyz . (5.5) 



^The renormalization of the coupling constants and masses can be made by just the standard way, 
by introducing the appropriate counterterms of renormalization in our model Lagrangian SB]. 
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For convenience we also choose parameters such that a = 1, G = 1 and A^. = 1, which 
means we consider rfi^ = 6|m^| and g"^ = X^f, = X^, which is a choice consistent with the 
considerations used in the previous Section, Eq. ( [4.4| ). We also take as base values for the (di- 
mensionless) dissipation coefficients the values {rj} = {TjxjTjyjTixy ) = (1/120, 1/240, 1/200)0, 
for which the dynamics displayed by ( |5.5| ) happen in the weakly damped regime. We also 
make the special consideration that the temperature of the system, which is kept fixed at 
a value T, is large enough such that the high temperature approximation involved in the 
Markov limit for the dissipation kernels is valid, but smaller than the critical tempera- 
ture, Tc^, for phase transition in the \& field direction, T^^ > T 3> F^, F^, m^, m^, where 
T^^c:.12\ml\/ig' + X^/2). 

We then numerically solve the dynamical system with initial conditions taken such that 
at t = the potential in the Lagrangian density ( p.l| ) is symmetrical in both field directions. 
The system is then evolved in time till the symmetry breaking in the \E'-field direction 
occurs. We then look for chaotic regimes as rj^, rj^ and rj^^p are changed. Our choice for 
initial conditions to numerically solving the (dissipative) dynamical system (|5.5| ) is as follows. 
At the initial time we consider 

(0„ 0„ V-c, 4) \t=o = (4$cr, 0, 0, 0) . (5.6) 

A typical result for the evolution of the fields in time is shown in Fig. 1, which already 
shows a highly chaotic dynamics prior to symmetry breaking. 

Following the method of Box-Counting |llO|-|T2i|l6| , around the initial condition ( p.6|) it 



is then considered a box in phase space (for the dimensionless variables) of size 10~^, inside 
which a large number of random points are taken (a total of 200.000 random points were 
used in each run). All initial conditions are then numerically evolved by using an eighth- 
order Runge-Kutta integration method and the fractal dimension is obtained by statistically 
studying the outcome of each initial condition at each run of the large set of points. Special 
care is taken to keep the statistical error in the results always below ~ 1%. The results 
obtained are shown in Table I for different values of dissipation coefficients. In this table we 



also show the uncertainty exponent e (see Ref. [|T^]), which gives a measure of how chaotic is 
the system. Qualitatively speaking we can say that the closest is e of zero, the more chaotic 
is the system. On the contrary, the closest is e of unit, the less chaotic is the system. We see 
clearly the effect of dissipation on the nonlinearities of the system. It can change fast from 
a chaotic to an integrable regime with a relatively small increase of the dissipation. Larger 
dissipations tend also to destroy fast the chaotic attractors. In Fig. 2 we show an example 
of the structure of the chaotic attractors in the ipc^Pc plane {y,w plane). 

VI. CONCLUSIONS 

In this paper we have studied the dynamical system made by the leading order effective 
equations of motion for a model of two coupled scalar fields. The chaotic behavior for 



^ This is consistent with the general expressions shown for the dissipation kernels Fi , F2 and F3 , 
which for A,^ = = g'^, satisfy the simple relation between then: 2{Fi + F2) = 5i^3. 
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(ensemble averaged) field trajectories has been demonstrated and we have also shown that 
in the overdamped regime chaos gets completely suppressed. Chaotic motion can only 
develop in the underdamped or very weakly damped regime, in which case enough energy 
can be exchanged fast enough from one field to the other, making both 0^ and ipc to fluctuate 
with large enough amplitudes, leading to a highly nonlinear behavior that then precludes 
the chaotic motion of the system. 

It seems also that we can indirectly associate the chaoticity in the system with the 
equilibration rates of the fields, in close analogy with the one found between the Lyapunov 
exponent and the thermalization rate in perturbative thermal gauge theory [T^. We note 
from the results obtained here and from the numerical simulations we performed, that the 
smaller is the fractal dimension (or the larger is e) the fastest the fields equilibrate to 
their asymptotic states by loosing their energies to radiation, which will then eventually 
thermalize. A clear assessment to this interesting point, however, must be carefully studied 
from the complete nonlocal EOMs, which then restores the time reversal invariance of the 
equations of motion. 

Finally we would like to point out that the kind of model we have studied here and its 
generalization to larger number of fields is natural to be found in extensions of the standard 
model with large scalar sectors. Physical implementations of the model can be found, for 
example, in particle physics or condensed matter models displaying multiple stages of phase 
transitions, in which case the dynamics we have studied here would be likely to manifest 
between any of that stages and, therefore, with consequences to the phenomenology of that 
models. The model studied here has also a strong motivation from inflationary models 
(hybrid inflation) |]18[. In special, in this context, a model Lagrangian of a similar form 



of the one we studied here has been studied in Ref. [^, showing the possibility of chaotic 
behavior during the final stages of inflation. However, the authors in |]TU[ make use of 
the classical equations of motion. It would be extremely interesting to assess the effect of 
quantum effects (and particle productions) and consequently dissipation also in that context, 
which is of fundamental relevance for the description of the process of reheating. In the 
particular study performed here, our results could apply instead to the description of the pre- 
inflationary stage, with possible contributions to the discussion of the fine-tuning problem 
of the initial field configuration in hybrid inflation |j20|. These problems are currently being 
examined by us and more results and details will be reported elsewhere. 
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FIGURES 




FIG. 1. An example of highly chaotic dynamics displayed by for two initial conditions 

given (in the dimensionless variables) by {x,z,y,w) = (4/\/6, — 10~^, 0, ±10"^). Time is scaled as 



r = 10y6|m2| t. 



-joU' '^.2 6 02 " 6'4 " '6'6 as ' 

y 

FIG. 2. The chaotic attractor for one reaUzation of 

the initial condition {x,z,y,w) = (4/\/6, — 10~^, 0, 10^^), for the first dissipation case shown in 
Table I. 
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TABLES 



TABLE I. The fractal dimension Df and the uncertainty exponent e for increasing dissipation 
coefficients (D = 4 is the phase space dimension). 

M = {Vx,Vy,Vxy) Df e = D-Df 

0.50 {r/} 3.96 0.04 

0.75 {r]} 3.90 0.10 

1.00 {r]} 3.43 0.57 

1.25 {ri} 3.23 0.77 

1.50 {;/} 3.10 0.90 
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